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, We present an analytical expression for the static many-body local field factor G+ (g) of a homo- 



o 



geneous two-dimensional electron gas, which reproduces Diffusion Monte Carlo data and embodies 
the exact asymptotic behaviors at both small and large wave number q. This allows us to also 
^ . provide a closed- form expression for the exchange and correlation kernel K,^c{r), which represents a 

^ ' key input for density functional studies of inhomogeneous systems. 

PACS number: 71.10.Ca, 71. 15. Mb 

The static charge-charge response function Xcio) of a paramagnetic electron gas (EG) can be written in terms of 
^ ' the Lindhard function Xo(<?) by means of the spin-symmetric many-body local field G+{q) through the relationship 
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Xc{q) = r, 7^T\^ — TT' 
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Thus G-i-{q) is a fundamental quantity for the determination of many properties of a general electron system. By 
I ] definition G-\-{q) is meant to represent the effects of the exchange and correlation hole surrounding each electrjOn in 
' O the fluid and is therefore a key input in the density functional theory (DFT) of the inhomogeneous electron gasH and 
C in studies oLquasiparticle properties (such as the effective mass and the effective Lande g-factor) in the electronic 
^ Fermi liquidH. 

For what concerns DFT calculations, a common approximation to the unknown exchange-correlation energy func- 
tional ii'xcM appeals to its second functional derivative, 
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where n is the average local density of the EG. The local field factor and the exchange-correlation kernel are simply 
related in Fourier transform by 

K^M = Jd\e-'^'^KUr)^-VgG+{q), (3) 

where d is the dimensionality of the system and Vg is the Fourier transform of the Coulomb potential e^/r. In what 
follows we shall only consider the case of two spatial dimensions, with d = 2 and Vq = 27re^/g. The corresponding 
three-dimensional case was discussed in Ref. Q . 
(~| ' A number of exact asymptotic properties of the static local field factor in two dimensions are readily proven. In 
O particular. 
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where kp = \/1'nn = y/2/rsaB is the Fermi wave number, rg — ^ nncLg is the usual EG density parameter with as 
the Bohr radius, kq = 7rrf/2 is the compressibility of the ideal gas in units of a^/Ryd, while k is the compressibility 
of the interacting system. By making use of the thermodynamic definition of k we can write 
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where ec(rs) is the correlation energy per particle. Once this function is known, it is possible to calculate A^. For 
the present purpose ec(f's) can be taken from the Monte Carlo data of JLef. 



The asymptotic behavior of G+ {q) at large q is also known exactly 



lim G+{q)^C+^ 



where C+ is proportional to the difference in kinetic energy between the interacting and the ideal gas, 

C+ = kp = 
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Moreover i?+ = 1 — g{Q), g{Q) being the value of the pair-correlation function at the origin. For g{0) we use the simple 
expression 
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1 + 1.372 rs + 0.0830 r2' 

which has been derived^ by an interpolation between the result of a low-r^ expansion, including the second order 
direct and exchange contributions to the energy in the paramagnetic state, and the result of a partial-wave phase-shift 
analysis near Wigner crystallizatiorL. rThis interpolation formula is in excellent agreement with many-body calculations 
based on the ladder approximationH-cl. 

In this work we fit the values of G+{q) originally obtained by Diffusion- Monte-Carlo (DMC) in Ref. in such a 
way as to obtain analytical expressions for both K^^{q) and K^^{r). Our formula for G^{q) reads 
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where q = q/kp and P+(g) is the polynomial P+{q) = <72 9^ + 94 q^ + 36 9^ + 58 9*- 

Some comments are needed in order to appreciate the correct physics which has been incorporated in Eq. (|1C 
(i) Our functional form embodies the exact asymptotic behaviors already introduced in Eq. and (^. (ii) The 
exponential factor e^"^-^^ ensures that G+(g) rapidly reaches the asymptotic behavior given by Eq. (^, a fact that is 
borne out by the DMC data at — iM- (iii) In the high-density limit (r-j 0) the term in square brackets tends to 
a two-dimensional Hubbard-like termllil, while the second and third terms tend to zero, (iv) The introduction of the 
high-degree polynomial P-i-(9) serves to reproduce the-jrich structure at intermediate wave number which is exhibited 
by G_|_ (q) as compared to the three-dimensional casellj. 

The only free parameters are contained in the last term in Eq. (^o|) and are fitted so as to minimize the differences 
from the DMC numerical results. For practical reasons it proves useful to have a continuous parametrization of the 
coefficients of the polynomial P+{q), which is at least valid in the range < < 10. We therefore propose the 
following: 
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1 + 0.8793 (r,/10)"-92i8 ' 
g2{rs) = 0.5824 (r./lOf - 0.4272 (r,/10) , 
g4{r,) = 0.2960 (r,/10) - 1.003 {rJlOf/'^ + 0.9466 {r,/lOf , 
ge{r,) = -0.0585 (r./lO)^ , 
g8(r.) = 0.0131 (r,/10)2. 

In Figure |l] we compare the fit given by Eq. (^M and ( pT| ) with the DMC data for = 1, 2, 5 and 10. In Figure ^ we 
show the local field factor G+ (q) as from Eq. (|10|) for various values of : the evolution from the low-r^ regime to 
the high-r^ one is clear. The fact that the highest peak in G+(q) occurs at = 5 is due to the behavior of C+{rs): 
this is a function that increases up to ~ 3.5, reaches a maximum and then decreases. Thus, the value of G+ at 
= 5 is larger than that at — 10. 

We turn next to the evaluation of K^^{r). ^From Eq. (||) and (10) the expression of the exchange-correlation kernel 
in real space (in Ryd) is readily obtained as 
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where Mi = -47r\/2C+/r5, M2 = ~2^/2B+/rs, M3 
—22 (;„ /vs- The function i^„(Q!, x) is given by 



-4\/2A+(l - e'^'/i°)/r,, M4 = %/2C+/r, and M. 
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where T{z) is Euler's Gamma function and iPi{a\b; z) is Kummer's function. In practice the function Fn{a,x) can 
be obtained via the recursive relation 
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where I„(z) is the modified Bessel function of order n. It is also useful to recall that dlQ{z)/dz — ii{z) and that 
dli{z)/dz — lo{z) — li{z)/z. 

In Figure |^ we show the exchange-correlation kernel K^^ (r) for various values of as from Eq. ( p^ (without the 
first term, which contains a two-dimensional S function). It is pleasing to note that no long-range oscillations are 
present in K^^{r). Notice that the M2 term in Eq. (12) diverges for kpr 0. 

In conclusion, we have presented an analytic parametrization of the local field factor entering the dielectric response 
of the two-dimensional electron gas in the paramagnetic state, incorporating the known asymptotic behaviors and 
giving an accurate description of the available quantum Monte Carlo data. We have obtained from it an analytic 
expression of the exchange-correlation kernel for density-functional calculations on inhomogeneous two-dimensional 
electronic systems. 
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